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SUMMARY 


The plane wave propagation, the stability and the rectangular duct mode 
problems of a compressible inviscid linearly sheared parallel, but otherwise 
homogeneous flow, are shown to be governed by Whittaker ’ s equation. The exact 
solutions for the perturbation quantities are essentially the Whittaker M-func- 
tions _j_^^^(4iTr] ) where the non-dimensional quantities t, n and 4 th have 

precise physical meanings. A niomber of known results are obtained as limiting 
cases of our exact solutions. For the compressible finite thickness shear layer 
it is shown that no resonances and no critical angles exist for all Mach numbers, 
frequencies and shear layer velocity profile slopes except in the singular case 
of the vortex sheet. 

INTRODUCTION AND BACKGROUND 

Studies on compressible free shear layers have not been extensive. In fact 
with the exception of the earlier work of Graham and Graham (ref. 1) who studied 
sound propagation through a finite linearly sheared layer in the low-frequency 
limit, it has been only recently that Blumen et al (ref. 2) obtained an exact 
solution for the stability of the shear layer with an hyperbolic tangent profile 
with the significant result that this shear layer is unstable to two dimensional 
disturbances for all Mach numbers whereas the vortex sheet is known to be un- 
stable only for M<2s/^ a result that cautions against modeling real sheared flows 
with vortex sheets, as has been the practice in a number of recent noise re- 
search studies, since, as the authors point out, even the long wavelength char- 
acteristics of finite thickness shear layers may be quite different from the 
corresponding properties of the analogous vortex sheet. In this study we con- 
sider sound propagation and stability in linearly sheared parallel compressible 
inviscid homogeneous flows. Work relating to the solutions of the pressure per- 
turbation equation has been that of Kiichemann (ref. 3) who also considered the 
stability of a boundary layer approximated by a linear velocity profile, the 
study of Pridmore-Brown (ref. 4) and that of Graham and Graham (ref. 1) . 

Kuchemann (ref. 3) obtained a formal series solution for the density per- 
turbation equation and he also arrived at a solution supposedly valid for large 
values of (our) parameter r\ = h His series solution, although it is given 

in a Clamber some and lengthy form, is correct but his asymptotic solution is in 
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serious error. Pridmore-Brown (ref. 4) solved the pressure perturbation equation 
in the short wavelength approximation, i.e., for large values of (our) parameter 

X = — K. His asymptotic solutions may also be in error. Graham and Graham (ref, 
4b 

1) studied the problem of a plane wave incident on a linear velocity profile free 
shear compressible inviscid layer. They used entirely a series solution of the 
density perturbation equation which they independently rederived apparently un- 
aware of the earlier work of Kuchemann. Because they used the series solution 
only, they were unable to give proofs for the range of the parameters t and n 
for which ordinary, total or amplified reflection occurs, although correctly 
identified the regions intuitively. More importantly and for the same reasons, 
they could neither prove the existence or non-existence of reasonances or crit^ 
ical angles, but for the case of "sufficiently thin — but not zero thickness — ■ 
shear layer" nor could they draw any conclusions for either the large Mach 
number M or large t cases. 

SYMBOLS 


a Speed of sound of the homogeneous fluid 

b Velocity profile slope of the shear layer 

c Disturbance phase speed (in general complex) 

Unit vector in the direction of the mean flow 
f,g Independent solutions of equation (1) 

k,kQ Wave numbers of incident wave 

m Second index of the Whittaker M-functions 

n Index of refraction 

n Wave normal unit vector 

Linear combinations of f and g in equation (5) 
q Parameter of transformation in equation (1) 

- 4 - ^ , 

w(r) z-component of the velocity perturbation 

x,y,z Rectanular coordinates 

Shear layer thickness 

A,B,C,D Functions of f, g and their derivatives in equation (7) 

A,B Functions of A,B,C and D in equation (9) 

H Heaviside function 

K Inverse of the x-component of the disturbance Mach number 

Disturbance vector Mach number 
d 

Mean flow vector Mach number 

M(z) Mean flow Mach number 

// // 

R.P. Real part of 
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^ l '^2 


T1.T2 

U(z) 


W 


a 

n 

0.00 

? 

X 

u i 
03 


Reflection coefficient 

In and out of phase component of R 

Transmission coefficient 

In and out of phase component of T 

Mean velocity 

Dependent variable in Whittaker Ts equation 

X- component of the wave vector of the pressure disturbance 

Non-dimensional variable in equation (2) 

Angle of incidence of plane wave 

Independent variable in Mhittaker's equation 
Non-dimensional parameter in equation (2) 

Perturbation velocity potentials 

Disturbance frequency 

SOLUTION OF THE PRESSURE PERTURBATION EQUATION 


In a homogeneous inviscid compressible parallel shear flow having a linear 
velocity profile in the z-direction only, i.e., U = U(z) = bz, it may be easily 
shown that, starting either from the linearized equations of motion or directly 
using the appropriate linearized form of the convective wave equation, the z- 
dependent part p(z) of the pressure perturbation p(?) is governed by the equa- 
tion : 

2 

where 

n = 7 - M, 4t = ^ K, and M = M(z) = ^ . (2) 

D a 

M is the local Mach number and K and w acquire the following meanings depending 
on the problem at hand: 

(i) Free Shear Layer : Propagation of a plane wave of wave vector k and fre- 

quency 0 ) impinging on the shear layer from a half-space (z<0) of relative 

rest and at an angle 0 measured from the z-axis (- — < 0 < + — ): 

2 — -- 2 

K = - — 

sxn0 

(ii) Free Shear Layer : Stability for assumed disturbances of the form 

p(r)= p(z) (a and c possibly complex): 

K = “ and 03 = ac 
c 
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(iii) Sound Propagation in Rectangular Ducts : Modes for assumed disturbances 
of the form p(r) = p(z) (k and o) real): 


ak 

2 

The transformation P “ W(^) with g — 4ix ^ reduces eguation (1) into 

Whittaker's equation for W, so that the two independent solutions f and g of 
equation (1) are: 


8 


= (4iT)"2~”M _ (4ixri^) 

it m 


m = - 


( 3 ) 


Where M are the \^ittaker M- functions and x,ri/K are defined following equ- 
ation (2) . It is easily shown that for x and n real^the functions f and g are 
also real with f being even and g odd functions of n whereas both f and g are 

even functions of x. Moreover series and asymptotic expansions for f and g are 

readily obtainable from the known properties of the Whittaker M’-functions (ref. 

5 and 6) . The series expansion agrees with the series obtained by Kiichemann 
(ref. 3) and Graham and Graham (ref. 1) although the form that results from equ- 
ation (3) above is not only more compact, but also faster converging. The asym- 
ptotic forms of f and g and their derivatives with respect to ti are obtained in 
terms of ^2m’fl with m = ± using Oliver's method (ref. 5 and 6) and they 

are in disagreement with both the results of Kuchemann (ref. 3) and Pridmore- 
Brown (ref. 4). This was expected as mentioned in the introduction since 

Kuchemann essentially seeking an expression for large n neglected 1 compared to 

in the last term of our equation (1) which is tantamount to setting x = O in 

Whittaker's equation, whereas Pridmore -Brown by applying hanger's method obtained 
only a non-uniform leading term of an asymptotic expansion in terms of Airy fun- 
ctions. These and other details may be found in reference 7. 


THE FINITE THICKNESS LAYER 


Plane Wave Propagation 

We consider the two-dimensional finite inviscid compressible shear layer of 
thickness with velocity profile 


U = bz^ z^z^ 

= bz 2 ;j^^z _>0 

= 0 0>z 


(4) 


in the (x,z) plane and a time-harmonic monochromatic plane wave incident from 
the z<0 half-space with wave vector it and wave momber k = in the (rest) frame 
of reference of the stationary fluid at z<0, in an otherwise homogeneous fluid 
in the entire (x,z) plane. The velocity potentials in the lower region (of re- 
lative rest) and the upper region of unifoim flow are: 


4 = 


R.P. 



±ik (xsin0+zcos0-at) , _ ±ik (xsin0-zcose-at) 
i +x\e 


, 2<0 
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= R.P.J+iTe 


±ik 


xsin0± (z-z^) I sin0 


1 


,z>z. 


U V 

where R.P, denotes "real part of"; the first term in equation (5) represents the 
incident wave coming from the half- space z^O and R and T are respectively the 
complex reflection and transmission coefficients for the velocity potential, and 
the upper signs are taken for ri>0-* In the middle (shear layer) region with the 
velocity profile equation (4) the pressure perturbation p(r) and the z-component 

of the velocity perturbation are given by: 


p (r) = p ^ ^ (n) sin 


w(r) = 


4xr) 


(n) cos 


k (xsin0“at) 

k (xsin0-at) 


+ P (n) cos 


( 2 ) . , . 
p^ (n)sin 


k (xsin0-at) 

k(xsin0“-at) 


where p^^^ (n) and p^^^ are linear combinations of the independent solution f and 
g, equation (3), of equation (1), i.e. 


(n) 


^ 11 ^ + ^ 12 ^' 


( 2 ) 


(n) = ^ 21 ^ + a 


"^ 22 ' 


(5) 


with a^j constants. Writing also R = R^ + the complex re- 

flection and transmission coefficients R and T respectively and applying the 
boundary conditions (continuity of the pressure perturbation p(r) and z-comp- 
onent w(r) of the velocity perturbation) at the interfaces z=0 and z=z^, and 

after separating real and imaginary parts in the resulting equations one obtains 
a system of eight linear algebraic equations for the determination of the eight 
unknowns a^ • R. , T. . After somewhat tedious but straightforward algebra the 

following expressions for the reflection and transmission coefficients are ob- 
tained: 


1 - r2 = ± Z2 

I sin20 I 


^ . ± ItaneiynJ - 1 

(A+B)^ + (CTD) 


(6) 


for |ri^|>l, - j <0£ R^, real, and 


- 1, T - - T„ - O, for |n-| |<1» - R- f T. complex 


2 - - 2 ' 


*This representation used by Miles (ref. 8) and Graham and Graham (ref. 1) is 
consistent with the radiation condition as postulated by Miles. Actually 
Sommerf eld's radiation condition does not apply for plane waves and the diffi- 
culties arising in such a case have been discussed by Lighthill (ref. 9) . At 
any rate these representations for the velocity potentials insure that the re- 
flected and transmitted waves are outgoing in a reference frame fixed in the 
upper fluid and are consistent with Miles's postulate and Ribner's intuitive 
picture (ref. 10) . 
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R^<1 

; Ordinary Reflection 



1 — 1 
II 
CM 

: Total Reflection 

(8) 


#>1 

: Amplified Reflection 



These regimes are quite analogous to the three regimes found by Miles (ref. 
8) and Ribner (ref. 10) for the vortex sheet case and intutiyely arrived at by 
Graham and Graham (ref. 1) . The above conditions in equation (8) imply that al- 
though the values of the reflection and transmission coefficients depend on the 
frequency u) and the velocity profile slope b, the dependence is through the 

angle of incidence 6, and the two parameters ^ sin0 and r)^= whereas 

the conditions for the three reflection regimes are independent or o) or b and de- 
pend only on tIt which is the Mach number —J: — of the x-component of the phase 

velocity of the incident wave 'front relative to the relative Mach number M^ of 
the two uniform flows confining the shear layer. 


The limiting case of the vortex sheet is easily obtained in the limit t^O 
(high-frequency or long acoustic wavelength limit) whereas the low-frequency or 
short wavelength limit is obtained by letting in equation (6); in the former 
case : 

1_r2 ^ +4 

which agrees with the results of Miles (ref, 8) and Ribner (ref. 10) , whereas in 
the short wavelength case 1-R^ becomes the Heaviside function H; 
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'f-Xx) 


i-r2 = h . 

Resonances and Critical Angles 

In this section we give a formal proof that in the amplified reflection re- 
gime, r)^<“l, there are no resonances and in the ordinary reflection regime, 

nj^>+l, there are no cricital angles. 

First it is easily deduced from equation (6) that excluding the singular 
cases of 0-K) or M-x» (in special ways) we may assume that neither A=G=0 or B=D=0, 
nor all four A,B,C and D may be zero simultaneously. It is next seen from equ- 
ation (6) that resonances exist if the denominator in the expression for is 
zero in the amplifying regime r]^<-lf i.e. , if 

A = - 0, A= A+B and B = C+D (9) 

with A,B,C, and D given by equations (7) . But A in equation (9) above is just 
the determinant of the coefficients a^^ in equations (5) of the system of the 

four equations determining . Since only the first equation of that system is 
inhomogeneous with right hand side proportional to x, a solution for the a^^j 
exists if and only if either A = A^ + B^ + 5 ^ O and t 0 or A = + B^ o'^and 

T^. It thus follows that equation (9) has solutions, only for x-KD, and this is 
precisely the limiting case of the vortex sheet; i.e. , resonances are possible 
only for the vortex sheet. One may also obtain the same result by algebraic 
manipulation of the general expressions for R^ and R^: 

2 2 2 2 

R B -A ^ -c — , R = 1 (ad - BC) , A = (A + B)^ + (C + D)^ (10) 

1 A ^ A 

for the case of amplified ref lection (lower signs) . 

For the critical angles we use the expressions in equation (10) with the 

upper signs (ordinary reflection n-, >1) • If critical angles exist, then R =R =0 

f f A ^ C 

and using equation (10) we may easily deduce that for B^^O, the ratio — = — 

B r) 

may then only attain the value -1 for zero reflection. But this implies that 
A+B=C+D=0 which is precisely the condition for the existence of resonances equ- 
ation (9) which we have just shown that do not exist for a finite thickness 
shear layer. 

THE FINITE THICKNESS LAYER 
Stability Considerations 

For the layer equation (4) the boundary value problem leads to the following 
equation: i(A + B) + (G - D) - O (11) 

where A,B,C, and D in general complex are given by equation (7) with K = f- 
with (sgn K) ommitted in the expressions for B and C. The roots of the above 
equation give the dependence of the phase speed c, or of the frequency ac, on 
the wavenumber a. For* temporal amplification a is real and positive, c and 
o)=acare real** For the neutral stability line however, in either case G=c^+ic^ 
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with c.=Or i.e., c is real, thus K is real and one may distinguish the cases 

X ^ ^ 

K^l, lr)|>^ correspoxidij^g "ho supersonic (upper signs) or subsonic (lower signs) 
disturbances and relative Mach numbers respectively. 

Comparing equation (11) and equation (9) we see that the resonances, were 
they to exist, would obey the system of equations A+B=0, C+D=0 whereas the 
neutral stability characteristiGs are determined by the system of equations 
A+B=0, C~D=0, with A,B,C,D real. Thus, in general one does not expect any con- 
nection between resonances and neutral stability eigenvalues except in the 
singular case of the compressible vortex sheet case which is discussed below. 

Special Case: The Compressible Vortex Sheet 


As before, we let x^O in equation (11) and (7) to obtain 

(1- + (i-Ti^)^ = o 


( 12 ) 


where x\ and K may be complex. Excluding the singular cases K->1 and ir^l# as well 
as K = -5^ we consider the case where the square root terms in equation (7) 
have the same signs, i.e., the cases K<l,|ri|>l or K>l,|ri|<l. The formal solution 
of equation (12) above is 


1 

K 


c 

a 





M± 

M^+4+4 Jl+M^ 

2 





which is in agreement with the results of Landau (ref. 11) . It is a matter of 
simple algebra to show that in order to satisfy the inequalities |ti|> 1 only 
the upper (minus) sign in the square root term in the above equation for the 
eigenvalues should be retained. Thus two neutral eigenvalues are not per- 
missible which is precisely the result of Miles (ref. 12) which he arrived at in 
a totally different way, namely by considering the vortex sheet stability as an 
initial value problem. It is finally worth noting that for the vortex sheet 
case the stability equation (11) for neutral eigenvalues, becomes C— D=0, whereas 
for the plane wave propagation case as we saw previously equation (9) for the 
resonances becomes C+D=0, since for the vortex sheet, t-h3, and A and B are 0(x) 
whereas C and D are 0(1) . Thus it is only for the vortex sheet that resonances 
and neutral stability eignevalues are given by the same equation i.e. C^ = D^. 
For the finite shear layer the roots of equation (9) and (11) the two equations 
are in general different. In fact, we have shown that although equation (9) 
may have real roots, there are resonances for the finite thickness compress- 
ible shear layer. 


Variable 


n = ± -M. 
K 


PHYSICAL MEANING OF r\ , T and 4Tn 

b 

1 c (o/k 


1/K is 


for plane wave propagation. 


sin0 a a 

stability and rectangular duct mode studies respectively, and we may thus write 
for T): -5- 

n = (M- - M ) -ef 
d r 
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Thus n is a relative Mach number measure i.e., it is the parallel to the 
mean flow component of the disturbance (phase speed based) vector Mach number 
M^, relative to the (relative ) mean flow Mach number and thus it is a measure 
of the components of the relative speeds of the disturbance and the mean flow 
in the direction of the mean flow. 

Parameter S. . acquires the simple meaning of a characteristic Strouhal 
b b 

number of the flow by writing: 

S = (Shear layer thickness) x iDisturbance frequency) ^ zo) ^ o) 

2 Mean Flow Speed bz b 

Parameters t and 4 tti : For propagation of a plane wave incident from a 

homogeneous half-space (z<o) at an angle 9^ with the z-axis, it is 

easy to show that the wave normals S are independent of x, i.e., that all wave 
normals of a given z-stratum are parallel. Thus defining an index of refrac- 
tion n = ( 1 +M^*n) = (1+Msin9).“ it is easy to show that n = l-MsinB^ = sinO^rif 

n = nsinOf). Using these relations we may write: 

^ sinS^rj^ = 4 tti, and 4 tti is the 

Vn|/„ I’ 

argument of the Whittaker M- functions in our general solutions f and g of the 
pressure perturbation equation. x itself also attains the simple meaning 

X = — K== Characteristic Strouhal Number 

b 1/K Parallel component of the disturbance Mach number 

CONCLUSIONS 


Local disturbance wavelength 
Relative refraction index change 


In this paper we have examined some aspects of plane wave propagation and 
stability of compressible inviscid homogeneous flows characterized by a linear 
velocity profile. The focus of this study has been on the search for exact 
solutions of the perturbation equations which bring forth the salient common 
features of all such parallel flows. The essential conclusions of this study 
are: (1) The z-dependent part of the pressure (or density) disturbance of 

such flows is governed by Whittaker * s equation with independent solutions 

(4iTn^) / with m 4 where fll are the Whittaker M-functions and ri/T, 

ix,±m 

and 4x11^ admit the following interpretations: 

n = (M^-M^) = Relative Mach niomber parallel to mean flow 


x= 


Strouhal number /^disturbance Mach number component 
/ in the direction of the mean flow 


4xri 


2 __ 


, V 

Vn /n 


Local disturbance v/avelength 


th j . 


^Relative refractive 
index change 


(2) Solutions to a number of other parallel flow problems may be obtained as 
limiting eases from our exact solutions. Such flows include the compressible 
vortex sheet (x-^O) , the incompressible vortex sheet (x-x», ir^O, (xn)-^0) y the in- 
compressible shear layer (x-><»,n^O, (xn) -finite) , and the short wavelength 
approximation of the compressible finite shear layer (x-^) . (3) The compress- 
ible finite thickness layer has no resonances and no critical angles for all 
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Mach n^IInbers, frequencies, shear layer thicknesses and shear profile slopes 
except for combinations of the singular values of O or «> for, w and b; two 
such combinations (b-x», z^-^0 but bz^ finite or oa-^) constitute the compress- 
ible vortex sheet case. 
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